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Abstract 



We analyze the density profiles and virial radii for a sample of 90 nearby 
clusters, using galaxies with available redshifts and positions. Each cluster 
has at least 20 redshifts measured within an Abell radius, and all the results 
come from galaxy sets of at least 20 members. 

Most of the density profiles of our clusters are well fitted by hydrostatic- 
isothermal-like profiles. The slopes we find for many cluster density profiles 
are consistent with the hypothesis that the galaxies are in equilibrium with 
the binding cluster potential. 

The virial radii correlate with the core radii at a very high significance 
level. The observed relationship between the two size estimates is in agree- 
ment with the theoretical one computed by using the median values of the 
density profile parameters fitted on our clusters. 

After correcting for incompleteness in our cluster sample, we provide the 
universal distributions functions of core and virial radii (obtained within half 
an Abell radius). 

Subject headings: galaxies: clusters of - galaxies: clustering 
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1 Introduction 



The distribution functions of observational cluster quantities, such as radii, 

velocity dispersions, masses, luminosities and X-ray temperatures, can pro- 
vide strong constraints both on cosmological scenarios and on the internal 
dynamics of these systems. Theoretical as well as observational estimates of 
these distribution functions are presently being debated. 

In previous papers (Girardi et al. 1993, hereafter GBGMM; Biviano 
et al. 1993), we have addressed the distributions of cluster velocity dis- 
persions and virial masses. Here, we deal with the spatial distribution of 
member galaxies in clusters. Existing relations between the linear sizes and 
other physical quantities of bound systems of galaxies have been investigated 
in the hterature, both observationally and theoretically (see, e.g., Oemler 
1974: Bahcall 1975; Peebles 1976; Dressier 1978; Sarazin 1980; Giuricin, 
Mardirossian, & Mezzetti 1984; West, Dekel, & Oemler 1987; West, Oem- 
ler, & Dekel 1989; Peebles, Daly, & Juskiewicz 1989; Lilje, & Lahav 1991; 
Cavaliere, Colafrancesco, & Scaramella 1991; Edge, & Stewart 1991; Inagaki, 
Itoh, & Saslaw 1992; Kashlinsky 1992; and references therein). 

Several size estimates are present in the literature: the core radius, the 
effective radius, the virial radius, the harmonic radius, the gravitational ra- 
dius, the mean radius, etc. (see, e.g., Bahcall 1977, and references therein). 
In the present paper we deal with virial and core radii. 

Relaxation processes naturally lead to the formation of physical cores in 
clusters of galaxies (see, e.g.. King 1966; Lynden-Bell 1967; Peebles 1970; 
Sarazin 1986). Generally, in order to fit the galaxy distribution of clusters, 
the following surface density profile has been assumed (see, e.g., Bahcall 1977, 
Sarazin 1980, and references therein): 

"^'^ = |i+(,vW 

where ctq is the central projected galaxy density and is the core radius. 
This model corresponds for r >> Rc to a volume — density p{r) oc r"^""-*-. 
Many authors (see, e.g., Bahcall 1975, Dressier 1978, Sarazin 1986) have 
used the King (1962) model to fit the cluster profiles; this corresponds to the 
model of Eq.(l) with a = 1. For instance, Bahcall (1975) found an average 
value of King — Rc — 0.12 ± 0.02 h^^ Mpc (the Hubble constant is Hq = 
h 100 km s~^ Mpc~^) on a sample of 15 rich regular clusters; on the contrary. 
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Dressier (1978) found an average value of King — Rc ^ 0.24±0.06 Mpc on 
another sample of 12 clusters. From the large compilation by Sarazin (1986), 
it may be noted that very different values of King Rc for the same cluster 
have been obtained by different authors. 

Core-radii have also been derived for the X-ray surface brightness distri- 
butions, by fitting these profiles to hydrostatic-isothermal models (sec, e.g., 
Abramopoulous & Ku 1983; Jones & Forman 1984; Lubin & Bahcall 1993; 
Bahcall & Lubin 1993; and references therein). An average value of ~ 0.05- 

0. 15 Mpc is suggested for X-ray core radii, although both smaller (Durret 
et al. 1993) and larger (Jones & Forman 1992) values have been found. 

Attempts to evaluate the core-radius of the cluster dark matter com- 
ponent have been made (see, e.g., Eyles et al. 1991; Gerbal et al. 1992; 
Kneib et al. 1993; Miralda-Escude 1992; Soucail 1992). The results are still 
preliminary. 

In this paper we use the hydrostatic-isothermal models in order to fit 
the cluster density profiles. In particular, this choice is corroborated by the 
recent results of Bahcall & Lubin (1993; see also Gerbal, Durret, Lachieze- 
Rey 1994), who claim that the so-called discrepancy" can be solved by 
fitting the galaxy density profiles with models having a — 0.7±0.1 in Eq.(l), 

1. e. p(r) oc r^^'"^"^*^'^ at large values of r, instead of using the classical King 
models. We remind that the /3— discrepancy is the difference between the 
observed values of Pspec, determined from cluster X-ray temperatures and 
velocity dispersions, and (3fit, determined from the gas and the galaxy density 
profiles in galaxy clusters (see, e.g., Sarazin 1986). 

The virial radius has been extensively used in the literature to evaluate the 
sizes and masses (via the Virial Theorem) of galaxy systems (see, e.g.. Limber 
& Mathews 1960; Smith 1980; Fernley & Bhavsar 1984; West, Oemler, & 
Dekel 1989; Perea, del Olmo, & Moles 1990; Pisani et al. 1992; Biviano et 
al. 1993). In the present paper we consider the 2D-projected virial radius 
(i.e. not de-projected by the classical value of 7r/2), 

^ ^^Tat (2) 

where N is the number of cluster members and Rij is the projected distance 
between the i and the j member. 

In this framework, we deemed it interesting to obtain estimates of cluster 
sizes using a large data-sample of 90 clusters, and we examine the relation- 
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ships between these sizes and other cluster properties. In § 2 we describe our 
data-sample; in § 3 we consider the method of analysis; in § 4 we provide the 
results and the relevant discussion; and in § 5 we list our conclusions. 

2 The Data-Sample 

We based our analysis on a sample of 90 nearby {z < 0.15) clusters, each one 
with at least 20 measured redshifts of cluster members within 1.5 Mpc . 
In this compilation we included some poor clusters (or rich groups) in order 
to enlarge our range of cluster-richness. 

The problem of foreground and background contamination may be severe 
when redshift information is not available for all the galaxies considered. In 
fact, different corrections have been adopted in the literature, leading to very 
different estimates for the sizes of the same clusters (see, e.g., Bahcall 1975; 
Dressier 1978; Fitchett & Merritt 1988; West, Dekel, & Oemler 1987). 

On the contrary, using observed redshifts and positions for all our cluster 
galaxies, we have been able to assign the cluster membership. Therefore, 
our sample should be considered quite free from contamination. The pro- 
cedure adopted to reject non-cluster members is similar to that described 
in GBGMM. However, in this paper we have decided to adopt a somewhat 
stricter selection criterion, which is described below. 

In order to reject interlopers in the clusters, we make use of the distance of 
a galaxy from the cluster center in the coordinates and in the velocity space. 
The center is defined as a biweight mean of the velocities, right ascensions 
and dechnations (see, e.g.. Beers, Flynn, & Gebhardt 1990, GBGMM). 

A first selection is performed to separate obvious groups which are pro- 
jected in the same cluster field. First, we identify all groups separated by a 
weighted gap larger than 4 in the velocity space. A weighted gap in the space 
of the ordered velocities is defined as the difference between two contiguous 
velocities, weighted according to their rank in the ordered distributions: the 
closer to the center of the distribution is the gap, the higher is its weight 
(for further details see, e.g.. Beers et al. 1990). We then select as our clus- 
ter the larger of these groups. Of the remaining galaxies, we subsequently 
reject all of them lying outside a circle of 3 Mpc from the cluster center. 
This rejection is made in an iterative way, recomputing the center after each 
iteration and thus redefining the 3 Mpc circle and the rejected galaxies. 



5 



until convergence is reached. 

This defines the main cluster body. Further selection requires the exclu- 
sion of galaxies more distant than 4000km s~^from the cluster mean velocity 
(as in Beers et al. 1991, GBGMM), and further rejection of galaxies sep- 
arated from the main cluster body by a weighted gap larger than 4 in the 
velocity space. The greatest difference with respect to GBGMM is in the 
last selection in the coordinate-space: we reject all galaxies more distant 
than 1.5 Mpc from the cluster center, once again in an iterative way. 
Since the size estimators are not as robust as the estimators of velocity dis- 
persions used in GBGMM, we believe that restricting the analysis to galaxies 
within 1.5 Mpc is a very simple but efficient method to reduce the prob- 
lem of interlopers, since the higher the aperture the larger the ratio between 
background/foreground galaxies and cluster members. Thus it is more con- 
servative to consider the inner cluster regions. 

Problems of magnitude incompleteness may cause errors in the size esti- 
mates (e.g.. Dressier 1978). In particular, if a deeper sampling of galaxies 
(in terms of magnitudes) has been made by observers in the central regions, 
the shape of the density profile will look steeper than it is in reafity. We 
have therefore analysed all our clusters with available galaxy magnitudes, by 
looking at the galaxy magnitudes as a function of clustercentric distance. We 
have noted that some clusters are sampled more deeply close to the center 
than they are outside. We have rejected the fainter galaxies of these clus- 
ters so as to erase any trend of the limiting magnitude with clustercentric 
distance. 

Yet another sort of bias may be present, if some regions in the cluster 
arc not covered obscrvationally (generally, in the outer parts). As for the 
previous case, the shape of the density profile will look artificially steeper 
than it is in reality. There is no easy way to eliminate this sort of bias a 
priori, but a signature of it may be the steepening of the density profile 
with increasing apertures. As a consequence, we look with more confidence 
on those results obtained from the central regions of our clusters, where 
problems of incompleteness should not be so strong (see §3). 

Yet another problem is the presence of substructures which may bias the 
determination of cluster sizes, and in particular of the core radii (see, e.g., 
Sarazin 1980). When the cluster has several density clumps, one density 
profile might not adequately fit the galaxy distribution; and large, spurious 
values of the core radius are expected when a single profile is forced to fit 



6 



all the data. In order to investigate the possible presence of subclustering in 
our sample, we performed the test of Dressier & Schectman (1988a). This 
test gives the probability, 1 — PdSi that some subclustering exists in the 
cluster (1000 simulations were run, according to the prescriptions of the 
above-mentioned Authors). Actually this method (like others available in 
the literature) is not efficient when the number of galaxies considered is low 
[N < 40). As a consequence, we considered the subsaniple of our clusters, 
with > 40 members and no evidence of subclustering, in order to check 
our results. As detailed in §4, our results do not seem to be biased by 
subclustering. 

The sample of the 90 clusters considered is listed in Table 1, together 
with the number A^ of galaxy members inside 1.5 Mpc from the cluster 
center, the Abell counts C, the richness R, the Bautz & Morgan and the 
Rood & Sastry classes (BM type and RS type), the probability P^s that 
observed structures in the cluster are chance fluctuations (i.e. low values of 
Pds indicate a high probability of subclustering) , and the relevant reference 
sources. 

3 The Method of Analysis 

In order to obtain the values of i?c, we fitted the observed galaxy distributions 
to hydrostatic-isothermal models on the hypothesis of spherical symmetry. 
Although clusters are not spherical symmetric (their flattening may indeed 
be quite large, see, e.g., Rhee, van Haarlem, & Katgert 1992), we prefer not to 
include ellipticity as a free parameter in our fits, because it would be poorly 
constrained in many clusters of our data-set, where we do not have enough 
data, and because inhomogeneous sampling during observing may result in 
artificial shapes for some clusters. We nonetheless performed Montecarlo 
realizations of distributions of 25, 50, 100 and 200 points extracted from 
King density profiles with varying Rc and ellipticity, and we analyzed how 
the derived values of Rc depend on the assumption of spherical symmetry. 
Not surprisingly, these values are intermediate between the values of the 
major and minor axis core-radii, and consistent with both within the formal 
errors. Therefore, we think that the assumption of spherical symmetry is a 
good working hypothesis. 

The fitting was performed using the Maximum Likelihood technique which 
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takes directly into account the positions of individual galaxies, and does not 
require any binning of the data (sec, e.g., Sarazin 1980). This method is 
very efficient, since it does not introduce any artificial scales, and it allows 
the simultaneous fitting of Rc and a. On the basis of the discussion of the 
previous section about interlopers, we assumed that the background galaxy 
density was negligible, and we set it at zero. 

We allowed Rc and a to vary from to 1 and from 0.5 to 1.5, respectively; 
we rejected values of Rc larger than 1 Mpc and values of a outside 
the above mentioned interval. Moreover, we did not consider size-estimates 
for clusters whose fitted profiles are rejected by a Kolmogorov-Smirnov test 
(hereafter KS test) at 90% of probability. However, fitting both Rc and a at 
the same time is not an easy task, especially when the data-set is poor. In a 
second stage of the investigation (see §4), we therefore decided to fix the value 
of a at the median value obtained in the two-parameter-fitting procedure. 
By reducing the free parameters to one, Rc, we were able to better constrain 
the size-estimates. 

The errors on our Rc were estimated by following the prescriptions of Avni 
(1976), while the errors on our Ry were obtained via the Jacknife technique 
(see, e.g., Efron 1979, Geller 1984). 

4 Core and Virial Radii 

We obtained both the virial radii and the best fits of cluster profiles with 
three increasing apertures, i.e. 0.5, 0.75, and 1 Mpc . 

In Table 2 we give the virial radii for our sample of clusters; in col.(l) we 
hst the cluster name; in cols. (2), (3), and (4) the values of -R^,o.5, Rv,o.75: 
Ry^i and their errors, computed within the increasing apertures 0.5, 0.75, 
1 Mpc for each cluster (if at least 20 members are contained within the 
respective aperture). 

In Table 3 we give the median values of the core radii i?c,o.5; Rc,o.75, Rc,i, 
and of the a's ckq.s, cko.ts cki for our clusters, computed within the apertures 
of 0.5, 0.75, and 1 Mpc , respectively. The numbers of cluster fittings, 
the core radii, and the values of a are listed in cols. (1), (2), and (3). The 
fittings were performed only if at least 20 members are contained within the 
respective aperture. Only the good fits were considered (see §3). 

The fitted profiles seem to be steeper when apertures larger than 0.50 
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h-^ Mpc are considered. This trend may be explained either by samphng 
incompleteness, or by a true difference between the profiles of the inner and 
the outer cluster regions (see, e.g. Bahcall & Lubin 1993). When we restrict 
our analysis to the inner regions, we note that the distributions of ao.s a-nd 
do not differ significantly, according to KS- and Sign-test; the same 
holds for /2c,o.5 ^-^d of i?c,o.75- The profiles seem to be well-defined within 
0.75 h^^ Mpc ; therefore, we restricted our analysis to this inner region. 

Fixing the a parameter in the fitting procedure makes it possible to reduce 
the formal errors on the derived values of the core radii. Our median value 
of ao.5 is O.SIqj. It does not differ, within the errors, from 0.7 ± 0.1, which 
is the a- value suggested by Bahcall & Lubin (1993) (see also Bahcall 1977) 
to resolve the /3— discrepancy for clusters of galaxies. So, we chose the value 
of 0.8 in our a-fixed models. 

In Table 4 we list the core radii obtained with the fixed a = 0.8 profiles; in 
col.(l) we list the cluster name; in col. (2) we list the values of Rc,o.5,a=o.8 ^-^d 
their errors, computed within the apertures of 0.5 h^^ Mpc , and in col. (3) 
the corresponding confidence levels (in percent) Po.5,a=o.8 of the profile fittings 
(obtained by using the KS-test); in col. (4) we hst the values of -Rc,o.75,a=o.8 
and their errors, computed within the apertures of 0.75 Mpc , and in col. 
(3) the corresponding confidence levels (in percent) Po.75,a=o.8 of the profile 
fittings (obtained by using the KS-test). 

As usual, the fits were performed only on samples with at least 20 galaxies 
contained within the aperture considered. 

The median values and the distributions of Rc,o.5,a=o.8 a-nd Rc,o.75,a=o.8 
do not differ, according to the KS- and Sign-test. This suggests that our 
estimates of the core-radii are quite stable within half an Abell radius. 

In order to detect a possible effect of subclustering on the radii estimates, 
we considered the subsample of our clusters with A^" > 40 members (for the 
reasons outlined in §2). On this sample we compared the distributions of Rc 
(and R^) for clusters with and without evidence of subclustering: there is no 
difference, according to the KS-test. We caution that this result cannot be 
considered neither general, since our cluster sample is not volume-complete, 
nor conclusive, since different tests for subclustering may yield different re- 
sults (see, e.g. West & Bothun 1990). In fact, a detailed analysis of 16 
well-sampled cluster (Escalera et al. 1994) has shown that the estimates of 
cluster virial radii do decrease when galaxies belonging to subclusters are 
identified and removed from the sample. However, the changes are not sta- 
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tistically significant in most cases, being important only when the cluster has 
a bimodal configuration, and this happens for only 4 out of the 16 clusters 
considered. In most cases the masses of the substructure are ~ 10 % of their 
parent cluster masses, so it is not surprising that subclustering does not have 
a very large influence on the size determination. 

The analysis of the errors of /?,. and R,. suggests that the dispersion of 
the values of Ry and Rc is partially intrinsic, and not mainly induced by the 
errors involved. 

4.1 Core Radii vs. Virial Radii, and Other Relations 

Both the core and the virial radius describe the galaxy distribution inside the 
cluster, so the existence of a relation between these two quantities is quite 
natural. This functional dependence can be easily derived when one knows 
the density profile. 

On the hypothesis that the models we use fit cluster density distributions, 
we obtain the following relation 



where A is the aperture considered. See the Appendix for details. 

In our data-sample the values obtained for the core and virial radii are 
very well correlated, at > 99.99 % significance level, for both the aperture of 
0.5 and 0.75 h^^ Mpc . Fig. 1 shows log -Rc,o.5,a=o.8 vs. log i?^, 0.5 (panel a), and 
log i?c,o.75,a=o.8 VS. log i?^,o.75 (panel b). This correlation is not induced by 
obvious observational biases (i.e., the radii do not correlate with the cluster 
distance, or with the number of cluster members we use, etc.). These relations 
remain practically unchanged if clusters with significant subclustering are 
removed. The theoretical curves derived from Eq.3 are superimposed on the 
data in Fig. 1. 

As can be seen, the observational relations between core and virial radii 
are in good agreement with the respective theoretical relations, described via 
eq.(3) and obtained with our {a = 0.8) profiles. This strongly supports the 
choice of these profiles in the description of cluster density distributions. 

Moreover, it is possible to compare the measured Ry with the corre- 
sponding values Rv,comp computed from the observed Rc and the respective 
apertures using Eq.3. correlates with Rv,comp at a significance level higher 




(3) 
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than 99.99 %. For the virial radii the regression bisector hne (see e.g. Isobe 
ct al, 1990) is 

Rv,0.75,comp = 0.02(±0.04) + 0.96(±0.06) • i?„,o.75- (4) 

Fig. 2 shows Rv,Q.75,comp versus Rv,o.75- Obviously, it is also possible to com- 
pare the measured Rc with the corresponding values Rc,comp computed from 
the observed and the respective apertures using Eq.3, and the result is 
quite similar. 

As a consequence of this analysis, since the estimate of R^ is usually 
more straightforward than that of Rc, one can evaluate Rc when i?„ and the 
corresponding aperture are known. 

We also investigated the existence of a correlation between (and Rc) 
and other optical cluster quantities, such as the Abell richness, the morpho- 
logical types by Bautz-Morgan and by Rood-Sastry, as reported in Table 1, 
and the robust velocity dispersions computed as in GBGMM: we found no 
strong evidence of correlation. 

The comparison between optical and X-ray core radii has been exten- 
sively discussed with controversial results (see, e.g.. Lea, Silk & Kellogg 1973; 
Smith, Mushotzky & Selermitos 1979, Abramopoulous & Ku 1983). Detailed 
optical analyses of few clusters (see, e.g., Ku et al. 1983, Henry & Henriksen 
1986) suggest that optical and X-ray core radii are comparable within obser- 
vational errors (see also David et al. 1990). We considered the X-ray core 
radii, R^, of Jones & Forman (1984); the median value of a in their fits is 
~ 0.7±0.1, which agrees with our optical a = 0.8. The distribution function 
of Rx does not differ (according to the KS-test) from that of our Rc,o.75,a=o.8, 
when the common subsample of 18 clusters is considered. Fig. 3 shows these 
distributions. 

4.2 Radii Distributions 

Our cluster set is not volume-complete. Therefore, in order to estimate the 
cosmic distributions for our cluster core and virial radii, we should normal- 
ize our observational distributions to a complete cluster sample. For this 
purpose, we considered the Cluster Catalogue of Abell, Corwin & Olowin 
(1989; hereafter ACO), which is nominally complete for clusters with red- 
shift z < 0.2 and richness class R> 1. 
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In GBGMM and Biviano et al. (1993) we used the Abell richness dis- 
tribution to normahze the distributions of cluster velocity dispersions and, 
respectively, of cluster virial masses. The normalization makes sense because 
the richnesses is available for all AGO clusters, and both velocity disper- 
sion and virial mass correlates with richness. Since there is no significant 
correlation between cluster size and richness, it is no obvious that a similar 
normalization would work in this case. On the other hand, the same lack of 
correlation suggests that the cluster size distribution should not depend crit- 
ically on completeness. We decided to perform the normalization to richness, 
anyway. 

The normalization was accomplished by randomly extracting 1000 values 
of Rc and Ry from the corresponding values measured in our clusters, in such 
a way as to reproduce the richness distribution of ACO's cluster sample (poor 
sample statistics forced us to consider clusters with > 3 as belonging to 
the same class). 

We considered the values of -Rc,o.75,a=o.85 listed in Table 4, having excluded 
the i? = clusters (AGO is complete for R> 1). The corresponding AGO- 
normalized cumulative distributions is shown in Fig.4; the median value of 
this distribution is 0.17lo!o5 Mpc (where the 99% interval hmits of the 
median are indicated). 

Since the virial radii depend on the apertures used, we cannot pro- 
vide a single cosmic distribution of virial radii. So, we chose the aperture 
0.75 Mpc (half an Abell radius) and the values of Rv,o.75 given in Table 2. 
The AGO-normalized cumulative distributions for Rvfl.75 are shown in Fig.5. 
The median value (with 99% interval hmits) is 0.671^:^6 ^"^ Mpc . 

5 Conclusions 

The cluster sample considered here is the most extensive one in the literature 
for the evaluation of virial radii and the fitting of galaxy density profiles, but 
it is not volume-complete. 

Most of our cluster profiles (~ 90%) are fitted by hydrostatic-isothermal- 
like profiles with a = 0.8; this result is, therefore, quite general. 

Our optical results are consistent with the hypothesis that many clusters, 
within the errors and the apertures considered, are satisfactory described by 
isothermal models. In fact, the slope of our best-fit profiles is similar to that 
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of X-ray cluster profiles (e.g., Jones & Forman 1984). Therefore we agree 
with the analysis of Bahcall and Lubin (1993), who suggested a value of 
a — 0.7 ± 0.1 in order to resolve the /3-discrepancy. 

This result may seem surprising in view of the mounting evidence that 
subclustering is ubiquitous in clusters (see, e.g.. Bird 1994). However, as 
discussed by West (1990), evidence for subclustering does not necessarily 
imply dynamical youth for the cluster. Classical tests for subclustering do 
not make any distinction between the case of simple geometrical projection of 
a group onto a cluster field, the case of a group that is falling into an otherwise 
relaxed cluster, and the case of real substructure, i.e. a true perturbation of 
the cluster dynamical potential. Moreover, subcluster masses are often an 
order of magnitude lower than their parent cluster masses (Escalera et al. 
1994), and therefore their influence on the main cluster properties need not 
to be huge. 

Our virial radii correlate with our core radii at a very high significance 
level. The relationship between these two size-estimates obeys the theoretical 
relationship which links the virial radius to the density profile with a = 0.8. 

We provide the distribution functions of i?„ and Rc obtained within a 
fixed aperture of half an Abell radius. The analysis of the radii errors allows 
us to claim that the dispersion of the values of the cluster radii is partially 
intrinsic. The median value of our core radii is 0.17 h^^ Mpc and the median 
value of our virial radii (within half an Abell radius) is 0.67 h^^ Mpc . 

We thank Craig Sarazin for kindly providing us with his FORTRAN code 
for the Maximum Likelihood analysis of density profiles, John Hill for giving 
us an electronic copy of his cluster data, and Michael West, the Referee, for 
his useful suggestions. This work was partially supported by the Ministero 
per I'Universita e per la Ricerca scientifica e tecnologica, and by the Consiglio 
Nazionale delle Ricerche (CNR-GNA). 
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Appendix 



Both our core and virial radii {Rc and respectively) describe the galaxy 

distribution inside the cluster, so the existence of a relation between these 
two quantities is quite natural. This functional dependence can be derived 
when one knows the density profile. 

The 3D virial radius -R„,3d of a distribution of N cluster galaxies of masses 
rrii on the sky plane, within a circular aperture of radius A, is given by the 
usual relation expressing the total potential energy of the system. Since only 
projected positions are known for the galaxies, a de-projection factor (usually 
7r/2) is used to deduce Rv^3d from its projected value Ry, which is derived 
from the relation 



miirij 




{Al) 

which becomes our Eq.(2) on the hypothesis that galaxy masses do not cor- 
relate with galaxy positions in clusters. This hypothesis is largely assumed 
in the literature and it is justified by the absence of pronounced luminosity 
segregation in galaxy clusters (see, e.g., Biviano et al. 1992). 

The relation (Al) describes the total potential energy of the 2D distribu- 
tion of masses. It is thus possible to evaluate R^ for a given surface-density 
profile a{r/Rc) {r < A), knowing the total potential energy W2d{A) of the 
disk of radius A. This energy is: 

W2d{A) = ^ 27r $(r) a{r/Rc) r dr; {A2) 
2 Jo 

in (A2) the factor 1/2 avoids counting twice the same couples of mass ele- 
ments in the disk, and the gravitational potential $(r) (using polar coordi- 
nates (r', 4>)) is: 

'•2^ (T{r'/Rc 
h Jo r'^ + r'2 — 2r r' cos ( 

The function I{r/Rc) = I{x) (where x = r/Rc) may be expressed, using 
the complete elliptical integral of first kind K{y) (see Binney and Tremaine, 
pg.73, 1987), as 



Hr) = -Gi i _^ d^dr'^-GR^I{r/R^). (A3) 



^ 4 /^/^^ ^ J, 
Jo X + u 



Au X 
{u + xY _ 



du, {A4:) 
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where u = r'/Rc, and K{y) is defined by 

K{y) = f ^ dt. {A5) 

From these relations it is possible to obtain the relation between and 
Rc, for a given profile and aperture: 

^ = ^^ rm.u W ^ y =RcF{A/R^). (AG). 
Jo i[x) a[x) X ax 

The function F{A/Rc) can be evaluated for any given surface density 
profile cr(r/i?c)- For the profile 

1 + {r/R,f\ 

F{A/Rc) is well represented, for 0.01 < {A/R^ < 100 and with an accuracy 
better than 3%, by the relation 

F(m, ^ 1.193 m 
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Captions to the Figures 
Fig. 1: 

log i?c,o.5,a=o.8 vs. log i?^,o.5 (panel a) , and log i?c,o.75,Q=o.8 VS. logi?^,o.75 (panel 
b). The superimposed curves are derived from Eq.3. 



Fig. 2: 

Rv,comp versus Ry. The superimposed regression line (Eq.4) does not signifi- 
cantly differ from the bisector line. 



Fig. 3: 

The cumulative distributions of -Rc,o.75,a=o.8 (solid line) and of the X-ray core 
radii of Jones & Forman (1984) for the common subsample of clusters. 



Fig. 4: 

The ACO-normahzed cumulative distribution of the core radii. 



Fig. 5: 

The ACO-normalized cumulative distribution of the virial radii, evaluated 
within half an Abell radius. 
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Table 1 (continued) 

Col.(l): Cluster name; 

Col. (2): Number of member galaxies within 1.5 h^^ Mpc from the cluster 





centre; 


Col. 


(3): 


Abell number counts, C; 


Col. 


(4): 


Richness class, R; 


Col. 


(5): 


Bautz-Morgan type; 


Col. 


(6): 


Roods-Sastry type; 


Col. 


(7): 


Probabihty of absence of substructure; 


Col. 


(8): 


Reference to velocities. 
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